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Abstract 

Stochastic lattice gases with degenerate rates, namely conservative particle systems 
where the exchange rates vanish for some configurations, have been introduced as simpli- 
fied models for glassy dynamics. We introduce two particular models and consider them 
in a finite volume of size £ in contact with particle reservoirs at the boundary. We prove 
that, as for non-degenerate rates, the inverse of the spectral gap and the logarithmic 
Sobolev constant grow as £ 2 . It is also shown how one can obtain, via a scaling limit from 
the logarithmic Sobolev inequality, the exponential decay of a macroscopic entropy associ- 
ated to a degenerate parabolic differential equation (porous media equation). We analyze 
finally the tagged particle displacement for the stationary process in infinite volume. In 
dimension larger than two we prove that, in the diffusive scaling limit, it converges to a 
Brownian motion with non-degenerate diffusion coefficient. 
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1. Introduction 

We analyze some models of stochastic lattice gases with hard core exclusion, i.e. systems 
of particles on a lattice A C Z d with the constraint that on each site there is at most 
one particle. A configuration is therefore defined by giving for each site x G A the 
occupation number r\ x G {0, 1}, which represents an empty or occupied site. The dynamics 
is given by a continuous time Markov process, which allows the exchange of the occupation 
numbers across a bond {x,y} of neighboring sites x and y with a rate c x ^(rf) depending 
on the configuration 77. The simplest case is the symmetric simple exclusion process 
(SEP), in which c x>y (rj) = 1. We instead consider processes in which there are some 
constraints in order for the exchange to be allowed, i.e. the rate c Xty (r)) degenerates for 
some configurations 77. We mention that Glauber-like, non-conservative, models with 
degenerate rates have been studied in [1]. We shall discuss conservative stochastic lattice 
gases with degenerate rates both in infinite volume, A = Z d , and on a bounded domain. 
In the latter case we shall however allow exchanges of particles with external reservoirs 
at the boundary. 

This kind of interacting particle systems, called kinetically constrained lattice gases, 
have been introduced in the physical literature as simplified models for some peculiar 
phenomena of the "glassy" dynamics [24]. Let us recall the physical problem [9]. Ex- 
perimentally, a glass can be obtained by cooling a liquid fast enough in order to avoid 
crystallization. Below the melting temperature the liquid enters a metastable phase in 
which the relaxation time is rather long and increase dramatically if the temperature is 
further lowered. When the time to reach equilibrium becomes longer than the experimen- 
tally accessible time scales, the liquid freezes in an amorphous solid phase, which is called 
glass. A complete theoretical explanation of this "glass transition" is still lacking [22]. 
The first question to be settled is whether the glass is a new state of matter or a long 
lived metastable state, i.e. whether the dramatic increase of the relaxation time is due 
to an underlying equilibrium transition or is a dynamical phenomenon. Since no static 
divergent correlation length is detected and the structural properties show a very small 
temperature dependence [9], it is possible that a purely dynamical transition takes place. 
Thus a major goal is to understand the mechanism inducing the dynamical arrest which 
prevents the relaxation of the system. This should be related to other phenomena such 
as the stretched exponential decay of the structure function for temperatures close to 
the glass transition and the aging phenomena for supercooled liquids quenched to lower 
temperatures [8,28]. 

In this context, several approaches have been proposed. We mention the random 
first order scenario [16,17] and the mode coupling theory [12]; from them the following 
picture in the mean field approximation arises [8]. There exists a finite temperature Td 
such that for T > the configuration space consists of a single ergodic component, 
while for T < Td it is broken into many disconnected ergodic components. Moreover the 
relaxation time diverges as T | Td- An open question is how to go beyond this mean 
field picture, in particular to establish whether this ergodic/non-ergodic transition is an 
effect of the mean field approximation or is related to the behavior of real glasses. The 
kinetically constrained particle systems, originally introduced to get some insight on the 
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physical mechanism for the dynamical arrest [18], have recently [4,29] also been used to 
investigate how the mean field scenario has to be modified for short range models. 

The basic idea behind these models is that the motion of a molecule is inhibited by 
geometrical constraints due to the presence of surrounding particles. In particular, a 
molecule can be caged by neighboring ones and the cage must be opened to allow its 
motion. It is thus possible that these local constraints might produce, for a finite value 
of the density, a cooperative behavior inducing the slowing down of the dynamics. The 
kinetically constrained lattice gases, are therefore defined by choosing exchange rates c x>y 
which encode this cage effect. Despite their simplicity and the discrete character, they 
might capture, at least at a mesoscopic level, some of the key dynamical ingredients of real 
glasses. Moreover, they have recently been used [25] to model granular systems, which 
is another class of systems displaying a glass-like dynamical arrest, known as jamming 
transition [13,20]. 

One of the most studied among such models is the so called Kob-Andersen model 
(KA) [18]. Let m = 1, . . . , 2d — 1, the KA model is defined by choosing c Xjy (rj) = 1 if 
at least m nearest neighbors of x different from y are empty and m nearest neighbors 
of x different from y are also empty, c x ^ y (rj) = otherwise. Note that for m = we 
would recover the symmetric simple exclusion. Let us consider this model for d > 2 in a 
finite cube A C Z d of size £; it satisfies the detailed balance w.r.t. the uniform measure 
on the hyperplanes with fixed total number of particles. However it is not ergodic, for 
instance if d = 2 and m = 1 two fully occupied consecutive rows do not move. By the 
general theory of finite state Markov chain it is however possible to decompose the state 
space into irreducible components. A natural question is the asymptotic behavior, as 
i — > oo, of the probability pe(p) of the "maximal irreducible component", here p is the 
density of particles on the hyperplane. The ergodic/non-ergodic transition at a density 
p c G [0, 1] would correspond to Poo(p) = 1 for p < p c and Poo(p) < 1 for p > p c . An 
alternative definition is obtained by looking at the model directly in infinite volume. 
We note the process satisfies detailed balance w.r.t. each Bernoulli product measure p p , 
p G [0, 1] but there are other invariant measures, for instance some concentrated on single 
configurations. Denoting by P t the semigroup associated to the process, by the spectral 
theorem, we have that 



if and only if zero is a simple eigenvalue of the generator; in such a case we say that the 
process is ergodic in L 2 (p p ). If the process is ergodic for p < p c and not ergodic for p > p c 
we would then say that the ergodic/non-ergodic transition occurs at p c . 

In the case d = 3 and m = 2, numerical simulations of the KA model [18] suggested 
that such a transition, in the sense of the former definition, takes place at p c ~ 0.881. 
However in [29] it is shown that for d = 2, m = 1 and d = 3, m = 1, 2 we have p c = 1 
whereas for d = 2, m = 2,3 and d = 3, m = 3, 4, 5 we have p c = 0. The finite size 
corrections are also discussed in [29]; in particular it shown that for p close to one the 
thermodynamic limit is achieved for sizes of order L(p) = exp{c(l — p)" 1 } for d = 2, 
m — 1 and d — 3, m — 1 whereas L(p) = exp{exp{c(l — p)" 1 }} for d — 3, m — 2. 
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It is also possible to consider the KA model in a finite volume of size i in contact 
with particle reservoirs at the boundary, see [4, 19] for numerical simulations. The total 
number of particles is not anymore conserved and, for m < d—1, the system is ergodic in 
the whole configuration space and reversible w.r.t. the Bernoulli measure whose density 
is fixed by the reservoirs. However, the dynamical arrest is not ruled out since, as I — > oo, 
the speed of convergence toward the unique invariant measure might exhibit a crossover 
as a function of the density. A preliminary question is then the asymptotic behavior of the 
relaxation time, which might be defined as the inverse of the spectral gap of the generator, 
for i — > oo. Note that for SEP the relaxation time grows as I 2 uniformly in the density. 

Another physically relevant issue is the asymptotic displacement x(t) of a tagged parti- 
cle. Indeed, for supercooled liquids near the glass transition, a decrease of the mean square 
displacement E(x(t) 2 ) is experimentally detected [9] with the possibility of a crossover 
from a diffusive, to a sub-diffusive behavior related to the dynamical arrest. Let us denote 
by Aseif = Aseif(p) the diffusion coefficient of a tagged particle. For the KA model with 
d = 3 and m = 2 the numerical simulation in [18] for p < p c ~ 0.881 suggested the power 
law behavior -D se if = (p c — p) a , with a ~ 0.3. This evidence is however probably due to 
finite size effects; in fact in [29] it is argued that for each p e [0, 1) the diffusion coefficient 
is strictly positive. 

The main results of this paper are sharp asymptotics on the relaxation time and the 
diffusive behavior of the tagged particle for some kinetically constrained lattice gases. We 
are not able to prove these results for the KA model, but they will be obtained for two 
simpler models defined by different geometrical constraints. The first model is defined 
by an exchange rate c XjV which vanishes if the two neighbors of the bond {x, y} along 
the direction y — x are both occupied. In the second model the exchange rate c XjV across 
the bond {x, y} vanishes if more than half of its neighboring sites, i.e. more than 2d — 1 
neighbors, are occupied. These models are in the same spirit of KA, but the degeneracy 
of the rates is not comparable. Indeed, there are exchanges allowed for KA and forbidden 
for our models as well as the converse. There is however an important simplifying feature 
of the models defined above, which plays an essential role in the rigorous analysis: it is 
possible to construct a finite cluster of empty sites which, uniformly in the configuration 
on its complement, can be shifted using only allowed exchanges. The results of this paper 
show that the models above defined behave essentially as the simple exclusion process and 
therefore they are not really appealing as models for the glass transition. The asymptotic 
as p | 1 of the relaxation time and of -D S eif(p) are however different from SEP. Of course, 
an interesting issue is whether the results obtained in this paper holds also for the KA 
model, in other words if the simplifying feature mentioned above is only a technical need 
or the behavior of KA is essentially different. As argued in [29], to which we refer for a 
further discussion on this point, it is expected that also KA is ergodic in /^(/ip) and the 
tagged particle diffusion is not degenerate for each p e [0, 1). The asymptotic of relaxation 
times as £ — > oo and the behavior of D sc \{(p) as p f 1 might however be different; their 
analysis appears to be a more difficult task. 

Another relevant issue is the macroscopic behavior of the kinetically constrained lattice 
gases. For non-degenerate rates, the hydrodynamical limit [27] states that if the initial 
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condition has a density profile then, under a diffusive rescaling, at later times we still 
have a density profile which can be obtained from the initial one by solving a parabolic 
equation. For SEP this is simply the heat equation. If the rates degenerate a natural 
candidate for the hydrodynamic limit is a parabolic equation of porous media type degen- 
erating when the density approaches one; an analogous result has indeed been proven for 
a model in which the occupation number rj x is a continuous variable [11]. For kinetically 
constrained lattice gases there is however a serious obstruction: the hydrodynamic limit 
cannot hold for any initial condition admitting a density profile. Consider for instance the 
models described above in d = 1 and take an initial configuration given by a sequence of 
two occupied sites and one empty at the left of the origin and three occupied sites and one 
empty at the right. This configuration is invariant for the microscopic dynamic, however 
the associated density profile evolves diffusively for the putative macroscopic evolution. 
Indeed, the initial density profile is 2/3 at the left and 3/4 at the right, it is thus bounded 
away from 1 and therefore is not affected by the degeneracy. On the other hand such phe- 
nomenon is somehow exceptional and we expect a hydrodynamic behavior for a suitable 
large class of initial conditions. 

Outline and summary of results. 

In Section 2 we define more precisely the models we analyze and introduce the basic 
notation. 

In Sections 3 and 4 we consider these models on finite volume of size £ with reservoirs 
allowing particle exchanges at the boundary. By our choice of the rates the processes are 
ergodic and reversible w.r.t. the product Bernoulli measure whose density is fixed by the 
reservoirs. We then discuss the rate of convergence to this unique invariant measure. 

In particular, in Section 3, we show that the spectral gap shrinks to zero as £~ 2 
when I — > oo. An analogous result for the simple exclusion process on the hyperplane 
with fixed number of particles has been obtained in [23, §8] and has been extended to 
Kawasaki dynamics, under suitable mixing conditions on the invariant measure, in [6,21]. 
Our proof is based on a comparison argument with the Glauber dynamics reversible w.r.t. 
the same Bernoulli measure and the construction of a suitable path which enables to move 
a particle from the boundary to any site by using only allowed exchanges. 

In Section 4 we show that the above idea can be used also to show that the logarithmic 
Sobolev constant, which controls the exponential decay of the entropy, grows as C 2 . An 
analogous result for the Kawasaki dynamics, under suitable mixing conditions on the 
invariant measure, has been obtained in [7,30]. We then prove that, via a scaling limit, 
one can obtain the exponential decay of a macroscopic entropy associated to the porous 
media equation. 

In Section 5 we consider these processes on infinite volume, in such a case they are 
reversible w.r.t. the Bernoulli measure p p and we show they are ergodic in L2(p p ) for 
any p G [0,1]. We do not discuss the rate of convergence to equilibrium in this context, 
which for the exclusion process is algebraic with diffusive exponent d/2, see [5,10,14]. 
We analyze the displacement x(t) of a tagged particle for the stationary process. We 
recall that for SEP in d = 1 we have [3] E(x(t) 2 ) P3 \fi while for d > 2 the displacement 
x(t) satisfies [15] a central limit theorem with strictly positive variance, i.e. 
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By the ergodicity of the process in L 2 (/i p ), one can repeat the arguments in [15,27] and 
show that, under a diffusive rescaling, x(t) converges to a Brownian motion with diffusion 
coefficient -D se if = D se u(p) given by a variational formula. Of course, if d — 1 we have that 
D se it = as for the simple exclusion. By following the strategy in [26,27] we finally prove 
that for d > 2 and each p G [0, 1) we have D sc \f > 0. 

2. Definition of the models 

The spatial structure is modeled by the <i-dimensional cubic lattice Z d in which we let e^, 
i = 1, . . . , d be the coordinate unit vectors. We denote by x, y, z the sites of Z d and by 
d(x, y) := \x — y\ the Euclidean distance between x and y. Given A, B C Z d we then let, 
as usual, d(A B) = inf {d(x, y) , x G A,y G £>}. 

For A C Z d , the configuration space in A is Qa '■= {0, l}" 4 . If A = Z d , we drop it 
from the notation, namely we let Q := VLjd. We can regard a configuration t] G as a 
map from A to {0, 1}; the value t] x G {0, 1} is interpreted as the number of particles of 
the configuration rj at the site x. If rj G Qa and 5 C A we denote by i]b the restriction 
of the configuration rj to f^- Let A, B C Z d be disjoint subsets, A n I? = 0; given 
?7 G Qa an d £ G we let 7/£ G ^aub be the configuration such that (t]£)a = f] and 
{vOb — For 7] G Q and x G Z, d we denote by the configuration rj shifted by x, 
namely ($xV)y := Vy-x, D G Given a probability measure \i and a random variable / 
we denote by fi(f) the expectation of / w.r.t. \x and by fi(f; f) := p (f — (f)) 2 its variance. 

For rj G Q we let T x i2/ ?7 = r] x ' y be the configuration obtained from rj by exchanging the 
number of particles in x and y, i.e. 



[ T) z if z^x,y 

Analogously, we let T x rj = rf be the configuration obtained from 77 by flipping the occu- 
pation number in x, i.e. 



We let also T x , resp. T XtV , act on functions / : Q — > R as T x f(rj) := f(T x rj), resp. 



T x , y f(ri) '■= f(Tx,yV)- We Anally introduce VJ := T x / - / and V x , y / := T^^J - /. 



The models we consider are defined as follows. Given a positive integer £, let A : = 
[1, £} d fl 7j d , A c = 7j d \ A. We consider the continuous time Markov process on the config- 
uration space with generator 




(2.1) 




(2.2) 




(2.3) 



where Ltmik describes the exchanges of particles in bulk; it is given by 




(2.4) 



{x,y}cA 
d(x,y)=l 
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namely the occupation numbers at the sites x, y G A are exchanged with rate c x>y . On 
the other hand Abound, describes the effect of the particle reservoirs at the boundary of A; 
it is given by 

Abound./ (v) = £ C *(v)^xf(v) (2-5) 

x£A,y£A 
d(x,y)=l 

For p G (0, 1) and x G A we choose the flip rates c x (r]) as 

c x {rf) := (1 - p)r) x + p(l - r) x ) (2.6) 

Therefore Lbound. is the generator of the following process. A particle on the interior 
boundary of A leaves the system with rate 1 — p (in fact at higher rate at the corners) 
while particles enter the system, if the landing site is empty, with rate p. We emphasize 
that in (2.3) the boundary part of the dynamics is slowed by a factor l/£, as we prove 
below this is the minimal choice to obtain that the relaxation time diverges as £ 2 . 

We let G Q be the configuration in which all the sites are empty. We shall discuss 
two specific choices of the exchange rates c Xjy (jj), rj G Q\. The first is 

c (l) / n - = / if (jl^x-ei + (V^A-)x+2e t =2, I = 1, . . . , d ^ ? s 

x,x+ei \ I) | J otherwise ' 

namely the exchange across the bond {x, x + e{\ is suppressed if the neighboring sites 
in the i direction are both occupied. Note that, since (t]0^)x = for x G" A, exchanges 
across the bonds {x, x + Ci} such that either x — G A c or x + 2ej G A c are not suppressed. 
The second is 



c 



xl(v) '■- { z:d({z},{x,y})=l ( 2 -8 




£ (r/0 A c) z > 2d - 1 



otherwise 



namely the exchange across the bond {x,y} is suppressed if more than one half of the 
neighboring sites are occupied. Note that for d — 1 we have c^ y = c x ]y. We shall denote 
by Lj[\ respectively Ly the generator (2.3) with c x chosen as in (2.6) and c XtV = ci^, 
respectively c x ^ y = c x 2 ) y . 

Let /xa,p be the Bernoulli measure on fl\ with density p, i.e. p/^,p is the product measure 
on f^A with marginal PA :P (r] x = 1) = p. It is easy to check that the generator k = 1,2, 
is self-adjoint in L 2 (/iA,p), equivalently the rates satisfy detailed balance w.r.t. p\ tP . 

We note that the bulk dynamics Lb u ik preserves the total number of particles in A, 
but - since the rates degenerate - it is not ergodic on all the hyperplanes of Qa with fixed 
total number of particles. For istance, if d — 1, all configurations t] in which the distance 
between all the empty sites is three or more do not evolve. On the other hand, thanks to 
-Abound., it is not difficult to show that the generator L * , k = 1,2, is irreducible, namely 
there is positive probability of going from any configuration to any other. By standard 
theory on finite state space Markov chain, irreducibility of implies the uniqueness of 
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the invariant measure and that is a simple eigenvalue of L A . In Section 3 we prove a 
lower bound on the spectral gap of in L 2 (pa, p ) showing that for each p G (0, 1) it 
shrinks to as £~ 2 . 

In order to discuss the diffusive behavior of the tagged particle we need to introduce 
also the infinite volume dynamics. The configuration space is then Q = {0, 1} Z , a function 
/ : Q — > K is called a local function if it depends only on finitely many t] x . The generator 
of the process acts on local functions as 

£ (k) f(v)= E '•^'/)V., ;/ /(,,) (2.9) 

{x,y}cZ d 
d(x,y)=l 

where c- k \ k = 1, 2 has been defined in (2. 7), (2. 8), where now we have A c = 0. Note that 
are translationally covariant in the sense that c^l yx+y+e {i}yri) = c x k ^ x+e {r]). Moreover, 
for each p G [0, 1] &ndk = 1, 2, the generator is self adjoint in L 2 (p p ), where p p is the 
Bernoulli measure in Q with density p. In Section 5 we prove that, for each p G [0, 1] and 
k — 1, 2, the generator & k > is ergodic in L 2 (p p ), namely that is a simple eigenvalue. 

We consider the process r] (t) generated by and condition that at time zero the 
origin is occupied; we then tag the particle at the origin and denote by x(t) its position 
at time t. The pair is then a Markov process on the state space {(rj,x) G 

Q x 7h d : rj x = l} with generator 

y&L d 

d(x,y)=l , 2 1Q . 

+ E c y %)[F( V y>*,x)-F( V ,x)] 

{y,z}GZ d \{x} 

Let Qq := {rj G f2 : r] = 1} and p p be the Bernoulli measure on Q with marginal 
Hp,o(Vx = I) = p, x E Zj d \ {0}. We shall consider the process (r/(t),x(t)) generated by 
A^ k ' with initial condition x(0) = and 77(0) distributed according to /i Pi o- in Section 5, 
for d > 2, we prove the invariance principle for the position of the tagged particle, namely 
that ex(e~ 2 t) converges in distribution, as e — >• 0, to a Brownian motion with strictly 
positive diffusion coefficient. 

3. Spectral gap 

The spectral gap of the Markov generator La is defined as 

gap(L A ) := inf spec (-L A \ l" 1 ) 

where 1 ± is the subspace of L 2 (pa,p) orthogonal to the constant functions. Since f2 A is 
finite, by irreducibility of L A , we trivially have gap(L A ) > 0, we next discuss its asymptotic 
behavior as i — > 00. 
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In the case of non-degenerate rates c x ^ y = 1, namely for the symmetric simple exclusion 
process, the aforementioned problem of non ergodicity on hyperplanes is not present. Let 
Vk,n(ji) '■= ^a,p(v\ XLeA^ = n ) be the canonical measure with n particles. In [23, §8] 
it is proven that, considering Lbuik with c x, y = 1 on L 2 (^A jn ), we have gap(Lbuik) x £~ 2 
uniformly in n, here ae x be means there exists a constant C > such that C~ 1 be < ae < 
Cbi for any I > 0. In the case c Xty = 1 it is not difficult to prove, and in fact it is a 
corollary of our analysis, that also for L A on L 2 (/iA, P ) we have gap(L A ) x £~ 2 uniformly 
in p. 

(k) 

Our first results is a lower bound on the spectral gap of L\ , k = 1, 2. Let us define 

(k) 

the Dirichlet form associated to L\ as 



1 1 (3-1) 

E A*A,p[c?2(V., w /) 2 ]+ I Y, /iA,p^(V,/) 2 ] 
{i,s}cA xeA,y^A ) 

d(x,y)=l d{x,y)=l 

Theorem 3.1. For any p G (0,1) and k = 1,2 there exists a constant C = C(d,p,k) 
such that for any i and for any function f : — > R we have 

VAAf-J)<Cl 2 S[ k) p (f) (3.2) 

Remark 1. Thanks to the variational characterization of the spectral gap, the bound 
(3.2) is equivalent to gap(L^) > C~ x l~ 2 . Moreover, by letting := exp{tL^} be the 
semigroup generated by L\\\ we also have that (3.2) is equivalent to 



\P t ik) f - VaJ\\l 2M < ex P { " 11/ " 



L 2(/iA,p) 



(k) 

for any function / on Q\. Finally, let E^, p be the distribution of the stationary process 

(k) 

generated by L A , i.e. we take p,^. p as the starting measure of the Markov chain. Then 
(3.2) is equivalent to 

^(/KO));/^))) <exp{-^} Wf-^fWl^ 
for any function / on Q\. 

Remark 2. By taking as test function f(rj) = XLeA^ — p) cos |f an d using c x k ) y < 1 
for any x, y and k — 1, 2, a simple computation shows that for each p e (0, 1) there exist 
a constant C = C(d,k,p) such that gap(L^) < C£~ 2 . Hence gap(L A % ' ) ) x £~ 2 as in the 
case of the simple exclusion. 

Remark 3. As discussed in the introduction, the correct dependence of the spectral gap 
on the density p has some interest. For simplicity we discuss it only in the case k — 1, 
namely for the rates chosen as in (2.7). It is a corollary of our analysis that the gap goes 
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to zero as p f 1 as a power law of exponent between 1 and 2. More precisely the following 
two inequalities hold. There exists a constant G\ = Ci(d) such that for any p £ (0, 1) and 
any I we have gap(L^) > (1 — p) 2 Ci/i 2 . For each integer £ > 5, there exists a constant 
C*2 = C 2 (d,£) such that gap(L^) < C 2 (l — p). The lower bound follows from the proof 
of Theorem 3.1; the upper bound is obtained easily by using as test function f{rj) = r\ x 
with x £ A such that d{x, A c ) > 3. 

Remark 4- As previously discussed, the process generated by L hu \ k is not irreducible on 
the hyperplanes with fixed number of particles, Qa^n '■— {v £ &a : YLxeA-Vx — A^}. In 
the one-dimensional case, d = 1 (recall that in such a case = c^), is however not 
difficult to check that Lbuik is irreducible on the set 

^a,tv := {v G ^a,n : 3 x, y £ A, x ^ y , d(x, y) < 2 such that rj x = 7] y = 0} 

and L bulk satisfies detailed balance w.r.t. the conditional measure i>a,n(') '■= ^a, p (-\^a,n)- 
A natural question is then the asymptotic behavior of the spectral gap of Lb u ik in L 2 {i , a,n)i 
a reasonable guess is that for each A^ < £—2 we still have gap(Ltmik) x ^~ 2 - This conjecture 
is supported by the fact that ii N = |A| — 2 = £ — 2, i.e. in the highest density case, we 
have a single pair of neighboring empty sites which performs a random walk. We are not 
able to prove the above conjecture in general, but only in the trivial situation in which 
A^ < |A|/3 = £/3. In such a case, for £ large enough, we have = Qa,n] therefore the 
statement follows easily by a comparison with the exclusion process with long exchanges, 
see [23, Lemma 8.1], and a minor modification of the argument in Lemma 3.2 below. 

The key step in the proof of Theorem 3.1 is the following lemma, which is in the same 
spirit of the path lemmata in [6,7,21,23,30]. 

Lemma 3.2. For k = 1, 2 and each x £ A let 



/iA.pMV,/) 2 ) < A[£ ^A,p(^+e 1 (V^ +e J) 2 )+ ^A c v(V y f) 2 )} (3.3) 



for any x £ A and any function f on Q\. 

Postponing the proof of the lemma above, let us first show how it implies, together 
with a comparison argument with Glauber dynamics, Theorem 3.1. 

Proof of Theorem 3.1. Let us introduce the product (Glauber) dynamics in defined 
by the generator 



U^ k) := {y £ A : 1 < Vl < x x - 1 , \ Vi - x-\ < 

where r^ := and r^ := 1. Then, for each k = 1,2 
constant A = A(d, k, p) such that for any £ 



rW,z = 2,...,ci} 



and p £ (0,1), there exists a 



yeul 



V- 3/1 = 1 




(3.4) 
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where c x has been defined in (2.6). The generator is self-adjoint in /^(/M.^); since it 
is a product dynamics, it is immediate to check its spectral gap is 1. For each function / 
on fl\ we thus get 



E ^p( c S + e 1 (v,, s/+e j) 2 )+ Yl ^(s(v,/) 2 )} 

<^(2r« + l)^ J] /iA,p(cS(V„/) 2 ) + ^p(^(V,/) 2 )} 

{x,j/}cA x£A,y£A 
d(x,y)=l d{x,y)=l 

(3.5) 

where we used the variational characterization of the spectral gap, Lemma 3.2 and ele- 
mentary inequalities. We thus get the bound (3.2) with C = A(2r^ + l) d_1 . □ 

We are left with the proof of the lemma. The basic idea is to use first Abound, to 
empty a few sites at the boundary. Then - via careful moves - we show how this cluster of 
holes can be shifted, using exchanges with non zero rate, and used to flip the occupation 
number in x. Finally we shift the cluster back to the boundary and use again Lb OU nd. to 
reconstruct the initial configuration near the boundary. 

Proof of Lemma 3.2. We discuss first the case of k = 1 which corresponds to the rates 
(2.7). We assume also that X\ > 4, otherwise the proof is much easier. 

Given rj E fl\ and x G A let us define S x rj G fl\ as the configuration given by 

if y = (l,x 2 , ...,x d ) 

if y= (3,x 2 ,...,x d ) v ' 




otherwise 



Moreover, for y G Z d we define 



± y -~ 1 y+e 1 ,y+2e 1 - L y,y+e 1 -L y-e 1 ,y /q n\ 

rpR rp rp rp ■ < J 



- V,y+ei 1 y+e 1 ,y+2e 1 1 y+2e 1 ,y+3ei 



Note that T y L moves the occupation number in y — e\ to y + 2e\ while the configuration 
in y,y + ei,y + 2e± is shifted by one in the direction —e\. Analogously moves the 
occupation number in y + 3ei to y while the configuration in y, y + e±, y + 2e\ is shifted 
by one in the direction e\. 

For x G A, we let 7^ := (iei,X2, ■ ■ ■ , Xa), i = 1, . . . X\ — 1 and define 

S xV ■= (7* • • ■T^_ 3 )T x _ eux (T^_ i ■ -TflStr, (3.8) 

It is not difficult to check that (S x rf) y = if y = 71 or y = 72, (S x ri)y = r] x if y = 73, 
(5*x?7)y = 1 - rj x if 2/ = a, and {S x rj) y = rj y otherwise. 
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For x G A with x\ > 4 we write 

V«/(»?) = [/fo* ) - /(F,i7)] + " /(Stf)] + [f(S,v) ~ /(i?)] (3.9) 

We start by considering the second term in decomposition above; we claim that 



XI — 1 



x; ^wp.w [/(^^) - /(^)] 2 < i8(i - P )- 2 (2x 1 - 7) x; ^,p[c« 7i+i (v 7ii7i+1 /)- 



To prove the above bound, let us introduce the path 



(3.10) 



i = 

il ) T x _ eux (T^ 4 ■ •■TflS.ri < = xx - 3 

K,.- " " " ^-3)^,4^1-4 -1%)S*V • = *i - 3 + 1, . . . ,2a* - 7 

(3.11) 

By (3.8) we have (2x l -7 = S x rj. Note also that, for each 77 e Qa and < i < xx — 4, 
the configuration Ci is guaranteed to be empty at the sites and 7«+2- Moreover, for 
each i| 6 Qa and < j < x± — 4 the configuration Czi-3+;?' is guaranteed to be empty 
at the sites ■y Xl -2-j and ■y Xl ~3-j- This will allow us to move the configuration Q to the 
configuration Ci+i using exchanges with non zero rate. 
By telescopic sums and Cauchy-Schwartz, we get 



(2xi-7 v 2 2xi-7 

E [/(Ci)-/(6-i)]j <(2x!-7) £ [/(Ci)-/(6-i)] 2 (3-12) 
i=i ' i=i 

We consider only the case 1 < z < xi — 4, the others are analogous. We then have 
Ci = T^Ci-ii recalling (3.7), again by telescopic sums and Cauchy-Schwartz, we get 



,12 

si— lj 



[f(Q-f (&-!)] = [/(T«Ci-i) - /(Ci- 

— 3^ [/(^7i>7i+ e i^7i+ e ii7i+2ei^7i+2ei.7i+3eiC«-l) ~~ /(2~7i+ei,7i+2ei^7i+2ei ,7i+3ei Ci 

+ [/(^7i+ei,7i+2ei^7i+2ei,7i+3eiCi-l) — /(^7i+2ei,7i+3eiCi-l)] 
+ [/(^7i+2ei,7i+3eiCi-l) - /(Ci— l)] } 

— 3^C 7ij74 _|_ ei (T 7i _)_ eii7i _ ( _2ei^7 i +2ei,7 i +3eiCi-l) [^7i ,7i+ e i / (^~7i+ e i ,7i+2ei^~7i +2ei ,7i+3ei Ci— 1 )] 

+ c 7i+eil7i +2ei(^7i+2ei,7i+3ei Ci-i) [v 

7i+ e ii7i+2ei i+2ei ,7i+3ei 

+ C 7i+2ei,7i+3ei(Ci-l) [V 7i+ 2ei,7i+3ei /(Ci-l)] } 



(3.13) 
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where we used that 



il,7i+3ei (Ci-l) c 7i +ei,7 1 +2ei (^7i+2ei,7i+3eiCi-l) 

= C 7i,7i+ei (^ 7 i+ e i:7i+2ei^ 7 i+2ei,7 i +3eiCi-l) = 1 

by construction of the path ^, see the remark below (3.11). 

In order to prove the bound (3.10) we consider only the last term on the r.h.s. of 
(3.13), the other can be analyzed in the same way. Given A C A, £ £ Q^, x & A, and 
g : f^A — > R we observe that, since fi^,p is a product measure and c x is given in (2.6), 



P^Av) VxC x {r]) g(VA\A0 = Yl VaAv) (1 - Vx) c x {rj) g{vA\ A 

= pO-- p)i*A, P ( L g\vA = £i) 

Recalling definitions (3.6) and (3.11), we thus get 



(3.14) 



Y^ A Av)c x (v) C 7i+2ei, 7 i+3ei(0-l) [ V 7i -f2ei ,74 +3ei 



i-1 



, 7 i+3ei 7i+2ei,7i+3ei 



(O-i) [v 7i+2ei,7i+3ei /(Ci- 

r?Gf2 

= P(l - P)|/ i A,p(4?+2ei, 7l +3 ei [V 7l+ 2 ei , 7i+ 3 ei /] 2 % = 1,7^ = ?7 7i+1 = ^ 7i+2 = o) 

+ ^A,p^C^+ 2ei , 7 i+3ei [V 7i +2ei,7i+3ei/] Vx = V^i = Vfi+i = 0; Vfi+2 = l) j 

(3.15) 

We now observe that for any positive function g : Ha ^ M we have 



7t+l 



0) 



PaA9) > PaAVx = 1, % = = % i+2 = 0) /iA, P (# I % = 1, P 7 , = *hi- 

+ PaAVx = Vn = Vji+i = 0, = 1) Pa, p (# I % = % = r? 7 ;+i = 0, 7] Ji+2 = 1) 
= P(! - p) 3 {pa, p (# I Vx = 1, »7 7i = % +1 = ^ 7l+2 = 0) 

+ Pa, p (# I Vx = = v^ +1 = 0, V^+2 = 1) } 

(3.16) 

so that from (3.15) we get 

5^ PAAv)Cx(v)^ + 2e 1 ,^+3e 1 (Ci^l) [V 7i+2ei,7i+3ei /(Ci-i 



< (1 P) 2 pA,p(c 7 .+ 2eij ^ + 3 ei 



[V 7< +2ei,7i+3ei/] J 



(3.17) 
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The bound (3.10) follows from (3.12), (3.13), and (3.17). Note that the extra factor 2 
comes from the return path, x\ — 4 < % < 2x\ — 7. 

Let us now consider the last term on the r.h.s. of (3.9); we claim that 
^2nA tP Cx(v) [f{S x rj) - f(r))] 2 

6 r 

71 ,72 L v 71 ,72 u 72 ,73 



- (1-p) 2 

(3.18) 

To prove the above bound let us define T+77 as the configuration given by (T£rj) x = 
1 and (T+rj) y = r\ y for y ^ x; analogously we let T~r] be the configuration given by 
(T~rj) x = and (T~7]) y = r\ y for y ^ x. Recalling (3.6) we then have 

f(S x rj) - f(rj) = rj x r/ 7l V 7l /(r/) + V 7li72 /(T~/7) + r/ 72 V 7l /(T 7li72 T"r/) 

+ ^72 , 73 / (^1-^71 ,72 -^71 ^7) ^~ ^71 ,72/(^72,73 ^71 ^1,72 -^71 ^7) 
"I" ^73 ^Jlf ( -^71 >72 -^72 ,73 -^71 -^71 ,72 ^71 V) 

+ (1 - r/ x .) (1 - 77 71 )V 71 /(?7) + V 7l!72 /(T+77) + r/ 72 V 7l /(T 7li72 T+r/) 

"I" ^72, 73/ (-^71 -^71172-^71 ? ?) "I" ^7i, 72/(^72, 73 -^71 -^71 ,72 -^71 7 ?) 
+ ^73 ^71 / (-^71 ,72 -^72 ,73 -^71 -^71 ,72 ^71 7 ?) 

(3.19) 

By using Schwartz inequality and the fact that the telescopic decomposition in (3.19) 
has been arranged so that all the exchanges have non zero rate, we get 



[/(^)-/(^)] 2 <6 ^%JV 7l /(77)] 



(l-^Xl-^JtV^r/)] 2 



+ r/ x r/ 72 [V 7l /(T 71i72 Tr/)] 



%)% 2 [V7 1 /(T 71;72 T+^)] 



+ % C 72,73 (^71 ^71 ,72-^71 ^) [^72,73/(^71 -^71,72 -^71 ^)] 

+ (1 - //^^(T-T^^T+r/)^^/^"^^^//)] 2 

+ Vx c -yi,j 2 (^72, 73 -^71^71,72 ^71 ? ?) [^7i ,72/ (-^72 ,73 -^71 -^7i ,72 -^71 ^)] 
+ ( 1 — 77 x ) C 7l ^ i72 (^72 , 73 ^ 71 T71 ,72 ^71 ^) [^7i ,72 / (-^72 ,73 -^71 -^71 ,72 ^71 ^ ] 
+ VxVyj, [^71/(^71,72^72,73^71-^71,72-^71^)] 

+ (1 — T] x )r]-y 3 [V 7l (/(T 71j72 T 72i73 X L^7i ,72^71^)] 



(3.20) 



By recalling the definition (2.6) we have 1 — r/ y < p 1 c y (r]) and 77^ < (1 — p) l c y {rf) 
for any 77 G f^A and y G A. We next estimate separately each term on the r.h.s. of (3.20). 
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(3.21) 



Let us consider only the last two terms, the others are easier. We have 

Pa,p(p) c x {v){ Vx P73 [^71 / (-^71 , 72 -^72 ,73 -^71 -^71, 72 -^71 PJ J 

+ (1 — T] x ) ?7 73 [V 7l (/(T 71i72 T 72j73 T^ i T 7i>72 T^ i ?7)] ^ 

< (1 - P) _1 {(1 - p) PA,p(^)pA, P (c 7l [V 7l /] 2 1 77^ = 1,?7 72 = 77 T3 = 0) 

+ P PA, P (1 - Px)pA, P (c 7l [V 7l /] 2 I ij x = r/ 72 = 0,r/ 73 = 1) 

< (l-p)- 2 /iA,p(c 7l [V 7l /] 2 ) 

where we used that, as in (3.16), for any positive function g : —> R 

fJ>A, P (g I Px = 1, P 72 = P73 = 0) + fJ,A, p (g I p* = ?7 72 = 0, ?7 73 = l) < ^ ^ ^ 2 Pa, p (p) 

By analogous computations for the other terms, (3.18) follows. 

Finally, to bound the first term on the r.h.s. of (3.9), it is enough to change variable 
77 1 — ► rf. Indeed, noting S x rf = S x rj, we get 

/WpMp) [f( V x ) - f(S x rj)] 2 = Y Pa,p(pK(p) [/(p) - f(S x r))] 2 (3.22) 

For x such that x\ > 4, the bound (3.3), with the constant A given by A = 180(1 — p) -2 
now follows from (3.9), (3.10), (3.18), and (3.22). The case in which 1 < x\ < 3 can be 
proven directly by the same steps leading to (3.18). 

In the case k = 2, namely for the choice (2.8) of the exchange rates, we give only a 
rough sketch of the proof, which is very similar to the case k = 1. Indeed, it is enough to 
define the configuration S x r], analogous to (3.6), as 

{0 if y e Qi((3,x 2 , • • ■ ,x d )) 

\-T) x if y = (l,x 2 , . . . ,x d ) (3.23) 
7] y otherwise 

where Q±(x) := {y G Z d : max i=1] _ )d \xi — yi\ < 1} is the cube of side 3 centered in x. 
The configuration S x r] can then be shifted using exchanges with non zero rates by 
means of a suitable path, depicted in Fig. 1 for d = 2. The proof is finally completed by 
the same arguments given for the choice (2.7). □ 



4. Log Sobolev inequality and entropy decrease for the porous media equation 

In this Section, we show how the techniques introduced in the proof of the spectral gap can 
be used to prove a logarithmic Sobolev inequality for exclusion processes with degenerate 
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Figure 1: The shifting path, o denotes sites guaranteed empty, • denotes l — r} x , ■ denotes 
an arbitrary occupation number, and — denotes the bond exchanged in the next move. 



rates. Moreover, by a scaling limit of the logarithmic Sobolev inequality, we deduce the 
exponential decrease of a suitable entropy for a nonlinear degenerate parabolic equation. 

We recall that a Markov process with reversible measure \i and Dirichlet form E is 
said to satisfy a logarithmic Sobolev inequality with constant c LS iff for any function / we 
have 

K /2ios ^)- cls£(/) (4i) 

It is well known, see e.g. [2], that (4.1) implies the hypercontractivity of the semigroup P t 
associated to the Markov process and the exponential decay of the entropy, namely, for 
any probability density / w.r.t. \x we have 

fjL (P t f log P t /) < exp { - ^ } n (f log /) 

(k) 

We first state the logarithmic Sobolev inequalities for the process generated by L A , 
k — 1,2. Recall the Dirichlet form E^ has been defined in (3.1). 

Theorem 4.1. For each k = 1,2 and p e (0,1) there exists a constant C = C(d,k,p) 
such that for any £ and for any function f : Q\ — > R we have 

"4 /2log ^)- c ^ (/) (4 - 2) 

Proof of Theorem J^.l. Let be the generator on introduced in (3.4); then, see 
e.g. [2], for any p G (0,1) it satisfies the logarithmic Sobolev inequality (4.1) with c LS 
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given by Ci(p) := (1 — 2p) 1 log[(l — p)/p] (we understand Ci(l/2) = 2) uniformly in £. 
By (3.3) and (3.5) the bound (4.2), with C = d{p) A (2r^ + l) 4 ' 1 (recall A = A{d,p, k) 
is the constant in Lemma 3.2), follows. □ 

Let us consider the following parabolic problem, called porous media equation, on 
B := [0, l] d with Dirichlet boundary conditions 

d t u(t, r) = V r ■ (D(u(t, r))V r n(t, r)) (t, r) G (0, oo) x B 
u(t,r)=p (t, r) G (0, oo) x dB (4.3) 

u(0, r) = (f(r) r G B 

where p G (0, 1), the initial datum tp G C(B; [0, 1]) satisfies (p(r) = p for r G <9.B, and the 
diffusion coefficient D(u) > is smooth and degenerates linearly for u — 1, namely the 
exists a constant 5 G (0, 1) such that 5(1 — u) < D(u) < 5 _1 (1 — u), u G [0, 1]. Since we 
assumed < ip < 1, by the maximum principle, we have that u G C(IR+ x B; [0, 1]). 

As discussed in the introduction, the equation (4.3) is the natural candidate for the 
hydrodynamic limit of the process with generator the diffusion coefficient D(u) 

would be given by a Green-Kubo formula [27, §11.2.2]. Note that the Dirichlet boundary 
condition is due to the particles' reservoirs. Although we do not prove any scaling limit of 
the microscopic dynamics to (4.3), we show how the logarithmic Sobolev inequality (4.2) 
implies the exponential decrease of a suitable entropy for the nonlinear evolution (4.3). 

Given p G (0, 1), we introduce the convex functional H p : C(B; [0, 1]) — > K + as 

HJu) := [ dr \u(r) log ^ + (1 - u(r)) log 1 ~ (4.4) 
Jb l P 1 - p J 

where we understand OlogO = 0. It is easy to show that H p is a Lyapunov functional for 
the evolution (4.3), moreover if u(t,r) is a smooth solution of (4.3) bounded away from 
and 1 we have 

-^-H p (u(t, -)) = fdr u(t, r) [1 - u(t, r)]D(u(t, r)) (v r log - u{t ^ ] ) * =: Q(u(t, ■)) 
dt J B V 1 — u(t, r)J 

(4.5) 

The following theorem, which states a "logarithmic Sobolev inequality" for the non- 
linear evolution (4.3) is easily obtained as a scaling limit of (4.2). 

Theorem 4.2. For each p G (0, 1) and 5 > there exists a constant C = C'(d, S, p) such 
that for any u G C 1 (B; [0, 1]) with u(r) = p for r G dB 

H p {u) < C'Q(u) (4.6) 

Remark. The inequality (4.6) can be proven directly by reducing it to the Poincare 
inequality for the Dirichlet Laplacian on B. The probabilistic proof given below, which 
somehow connects the evolution (4.3) to the microscopic process, shows additionally that 
the Lyapunov functional H p is the macroscopic limit of a relative entropy. 
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Proof. We shall prove the bound (4.6) for D(u) = 1 — u 2 ; the generic case of D de- 
generating linearly for u \ 1 follows by our hypothesis on the diffusion coefficient D. By 
truncation, it is also enough to prove (4.6) when u is a smooth function bounded away 
from and 1. 

We set e := t~ l and apply inequality (4.2) for k = 1 choosing f 2 = g £ where 

Note that ^% pu {t}) := ^A,p(v)9e(v) is a product probability measure on Qa with density 
profile w, namely ^% pu {f]x) — u(ex). By elementary computations which we omit, we have 
that the normalized relative entropy of ^ w.r.t. pi^ p converges to H p {u) as e — > 0, 
namely 

limeVpf&log gf ) = # p ( w ) (4.8) 



Moreover it is straightforward to check that 



hm£ d - 2 4%^) = QH (4.9) 



Let C(d,l,p) be the constant such that (4.2) holds for k = 1. By (4.8) and (4.9) the 
bound (4.6), with C = C(d, 1, p), now follows from Theorem 4.1. □ 

The exponential decrease of the "entropy" H p along the flow of the porous media equa- 
tion (4.3) follows from (4.5), Theorem 4.2, and a straightforward truncation argument. 

Corollary 4.3. Let u G C(R + x B; [0, 1]) be the solution of (4.3) and C = C'(d, 5, p) be 
the constant in (4-6). For each p G (0, 1) we have 

#pHV)) < e-' /C '# P (^) (4.10) 

for any t G M + and any tp G C(B; [0, 1]) such that (p(r) = p for r G dB. 

We have assumed that the diffusion coefficient D(u) degenerates linearly for u — 1. 
One can also obtain the exponential decrease of the entropy H p if D(u) x (1 — u) n , n 
a positive integer. This can be shown by introducing a microscopic model in which the 
exchange rate c XjX+ei (rj) is zero iff there exists j = 1, . . . , n such that f] x _j ei = rj x+ ^ + iy- = 
1. It is in fact possible to prove that the logarithmic Sobolev constant for such a model 
is of the order £ 2 . 



5. Diffusion of the tagged particle 

In this Section we consider stochastic lattice gases with degenerate rates in infinite volume. 
We first prove that, for each p G [0,1], the generator C^ k \ k = 1,2 defined in (2.9) is 
ergodic in Z^Ov); recall that \x p is the Bernoulli measure on Q. 

Proposition 5.1. For each p G [0, 1] and k = 1,2 we have that zero is a simple eigenvalue 
of the generator considered on L 2 (fi p ). 
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Proof. Let 

^\f) = \ £ ^Ki(^«/) a ] (5-1) 

{x,y}cZ d 
d(x,y)=l 

be the Dirichlet form of the generator C^ k \ To show zero is a simple eigenvalue of 
we check that Sp k \f) = implies / constant /i p -a.s. This is trivially true for p = 0, 1. 
For p G (0, 1), by De Finetti's theorem, it is enough to show that £p (f) = implies 
HpC^x,yf) 2 = for each {x, y} C Z d with d(x, y) = 1. 

We discuss in some detail the case k = 1. Let x G Z d and consider the bond {x, x + ej}, 
i = 1, . . . , d. For n — 1, 2, . . . we introduce the events 

#M := e : %+„ ei = % +(n+ i )ei = 0} B Xji := (J 

n>l 

By noting that p p (B x ^) = 1, we have 

oc 

fi p (V x , x+e J) 2 = p p {[V XtX+e J} 2 l B J < £/i P ([V x , x+ei /] 2 l B « J (5.2) 

n=l 

Let 7h := x + hei, h = 0, 1, ... ; given 77 G B xi we can find a path 77 = £ , . . . , (n = 

r] x ' x+ei where (j+i = (J h ' lh+1 for some h = 0,1,... and c$ jlh+1 ((j) = 1. It is in fact 
possible to construct a path analogous to the one introduced in the proof of Lemma 3.2; 
note that the two sites 7„, and 7 n +i are empty by the definition of the event B xi . Since 

4 1} (/) = implies 

(41^+1 i V ih,7h+i ff) = for any /i = 0, 1, . . . 

by telescopic sums and Cauchy-Schwartz in (5.2) we get p p (V XtX+ei f) 2 = 0. 

Recalling Figure 1 it is straightforward to modify the argument given above to cover 
also the case k = 2. □ 

We next discuss the diffusive behavior of a tagged particle. More precisely, we consider 
the process (77 (i), x(t)) with generator given in (2.10), initial condition x(0) = and 77(0) 
distributed according to p Pt0 , the Bernoulli measure on Q = {7/ G Q : 7/ = 1}. Let 
£(rj) := 'd- x (t)ri{t) be the process as seen from the tagged particle, we have that £(rj) is 
itself a Markov process on the configuration space VLq with generator given by 

4 fc) /(o= E c S(0(i-u[/(^°' y )-/(o]+ E c S(ov^/(o (5.3) 

j/GZ d {x,y}cZ d \{0} 
d(0,y)=l d{x,y)=l 

A straightforward computation shows that is self-adjoint in L 2 (p P fi)', moreover, 
by the same argument as in Proposition 5.1, it is also ergodic in Z^G^o)- We can therefore 
apply the same proof as the one given in [15,27] for non-degenerate rates. We get that the 
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rescaled position of the tagged particle, ex(e~H), converges in distribution, as e — > 0, to a 
d-dimensional Brownian motion with diffusion matrix 2 D^} { . Furthermore the diffusion 
matrix = D^ { (p) is given by the variational formula 

d(0,2/)=l 

+ E '30T.,,f((); 2 } 

Ry}CZ d \{0} 
d(x,y)=l 

(5.4) 

where r G I* 1 and ■ is the inner product in M. d . 

The main result of this Section is that, for d > 2 and each p G [0, 1), the diffusion 
matrix D^ { (p) is strictly positive as in the case of simple exclusion [15,27]. 

Theorem 5.2. For each d > 2, k = 1, 2 and p G [0, 1) i/iere exists a real c = c(d, k, p) > 
such that r ■ D^ f (p)r > cr ■ r for any r G M. d . 

Remark. As discussed in the introduction, the behavior of D^ f (p) as p f 1 has some 
interest. Note that for SEP it vanishes linearly. In the case k — 1, Theorem 5.2 will be 
proven with c = Cq(1 — p) u where Cq does not depend on p. However, by the same strategy 
and some extra efforts, it is possible to improve the lower bound to c = c' (l — p) 4 . An 
upper bound of the form D^ f (p) < Cq(1 — p) 2 l is easily obtained by using a constant 
test function / in (5.4). 

Let us fix a direction in M. d , say e\ and define the following subsets of Z d \ {0} 



R 



(i) 



{x G 7L d \ {0} : max |xj| = 1 , Xj = , % — 3, . . . , d} 



■= [ x G Z d \ {0} : xi = ±2, |x 2 | < 1, X; = 0,i = 3,. . . , d} 



and 



^ : = {x G Z d \ {0} : xi = 0, maxjxil < 3} 
: = {x G Z d \{0} : xi = ±1, maxjxil < 3} 



(5.5) 



(5.6) 



Given £ G f2o, we next define £ + ' ,( - fc ) as the configuration obtained from £ by exchanging 
the occupation numbers in R+l with the corresponding ones in R { X namely 

^ } * 1 ^ 4ei if x G itS ( } 



and 

if x G 
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We finally introduce the events 

Bf ■= {£ G n : ^ = , = 0} , Z3 (fc) := tff U £ W (5.9) 

and note that f G 8^ iff G #i fe) . 

Lemma 5.3. For each d > 2, k = 1, 2, and p G [0, 1) there exists a real a = a(d, k, p) > 
such that for any r G M d we have 

r ■ D^(p)r > (r • e x f \ J p pfi (d^) { [/(£+•-.«) - /(0] 2 



y=±l 



(5.10) 



Proof. We discuss in some detail the case k — 1. We note that if £ G Z^ 1 -* and y = ±1 
we have 

tiLm ~ ZyeA > Me C1) =0}(£) 
since /ip^iS'- 1 )) > (1 — p) 11 , by the same argument as in (3.16), from (5.4) we then get 

r ■ D<»(p)r > (r ■ e,f \ M { I „ ft0 (dfl £ cg> (0(V B , y /(0] 2 

{*,y}CZ d \{0} 
d(x,j/)=l 

+(i-p) u L p m\b {1) ) E % R , 1 )=o}(0[z/ + /(^e 1 ^ ei )-/(o] 2 } 

?/=±l 

(5.11) 

Let Ti, . . . ,T 16 be the chain of exchanges depicted in Figure 2, Tj exchanges the oc- 
cupation numbers in the bond 6j. Note that if £ G £>+^ the path (q := £, := Ti(q, 



Ci + 6 := 7i 6 d + 5 = is such that ^(C+J = 1, z = 1, . . . , 16. For £ G 



we 



define analogously Co := C-i Ci := ^i6Co > ■ ■ ■ > Ci6 := ^iCis = £ + ' which is such that 

<£_ 4 (Ci) = M = !,-••, 16- 
We then have 

[f(t'-' {1) )-m)]\ m (o 

< * B w(0[m + '-' (1) ) - m)] 2 + V)(o[/(e + '-' (1) ) - m)] 2 

16 16 

= v>(o[E/(c + ) - /(c + -i)] 2 + v^)E/(cr) - /(d)] 2 ( 5 - 12 ) 



i=l i=l 

16 



16 E {^(Ci) [v 6 j(cti)] 2 + 4!L(Ci) [v bl7 „j(Ci)] 2 } 



< 

j=i 
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Figure 2: Chain of exchanges Ti, . . . , Tie for k — 1. The picture represents sites in the 
plane e%, e2, • denotes site and — in the i-th figure denotes the bond exchanged by Tj. 
If the sites denoted by o are empty then the starting configuration is in B+ . In such a 
case 1,2, and 3 denote the occupation numbers in which, step by step, are moved to 
R+ by using only allowed exchanges. 

By integrating w.r.t. /i Pi o the above inequality and taking into account that in the chain 
of exchanges Ti, i = 1, ... ,16 each bond is used at most twice we get 

^o(B«) /^,o(^|B«) [/(e+.-.«)-/(0] a 

<64L p ,oK) J2 cgi(0[V,, w /(0] a 

{x,y}cZ d \{0} 
d(x,y)=l 

which inserted in (5.11) concludes the proof with a(d, 1, p) = 2~ 6 (1 — p) 11 . 

The case k = 2 is proven by the same arguments; in this case, for al = 2, the required 
chain of exchanges T±, . . . , T 35 is depicted in Figure 3. □ 

Proof of Theorem 5.2 (sketch). Thanks to Lemma 5.3, it is enough to prove that the 
right hand side of (5.10) is strictly positive for r-e\ ^ 0. By the variational formula (5.4), 
it can be interpreted as the self diffusion coefficient of a one dimensional auxiliary process 
which we next describe in the fixed frame of reference. 

The configuration space is {(y,Tj) G Z x O : ^~y ei V € S^}. Let y(t) E Z be the 
position of the tagged particle and r)(t) be the particles configuration. At time t = 



(5.13) 
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Figure 3: 


Chain 


of exchanges T 1; . 





figure denotes the bond exchanged by Tj. 



2. • denotes site and — in the 2-th 



the tagged particle is at the origin, y(0) = 0, and 77(0) G B^ is distributed according 
to /i p ,o( • |£>( fc )). Then the tagged particle jumps to the right, resp. left, with rate one if 

$-y(t)eiV(t) £ B+\ resp. if fi_ y (t) ei ri(t) G B_. Moreover, with rate one, r](t) is exchanged 

to $ y (t) ei [($-y{t)e 1 v(t)) + ~'^] > namely the occupation numbers in fiy^^R^ are exchanged 

(k) 

with the ones in l d y (t) ei R+ ■ 

The proof of the Theorem can now be completed as in the case of non-degenerate 
rates, see [27, II.6.3], by showing that there exists a real c > such that for any t > and 
rj G B^ we have E( 0iT) ) (y(t) 2 ) > ct. Here E( 0)f? ) denotes the distribution of the auxiliary 
process with initial condition (0,77). □ 
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